The Arp 220 merger on kpc scales 
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ABSTRACT 

For the first time we study the Eastern nucleus in greater detail and search for the more 
extended emission in the molecular gas in different CO line transitions of the famous ULIRG 
Arp 220. Furthermore we present a model of the merger in Arp 220 on large scales with the 
help of the CO data and an optical and near-infrared composite HST image of the prototypical 
ULIRG. Using the Plateau de Bure Interferometer (PdBI) we obtained CO (2—1) and (1—0) 
data at wavelengths of 1 and 3 mm in 1994, 1996, 1997 and 2006 at different beam sizes and 
spatial resolutions. The simulations of the merger in Arp 220 were performed with the Identikit 
modeling tool. The model parameters that describe the galaxy merger best give a mass ratio 
of 1:2 and result in a merger of ^6x10^ yrs. The low resolution CO (1—0) PdBI observations 
suggest that there are indications for emission ~ lO^'towards the south, as well as to the north 
and to the west of the two nuclei. 

Subject headings: Galaxies: interactions - Galaxies: kinematics and dynamics - Galaxies: nuclei - Radio 
lines: galaxies 



1. Introduction 

With an infrared luminosity of L = 1.4 x 10^^ Lq 
at a distance of Dl = 78 Mpc {Vt\ = 0.7, VLm = 0.3, 



and Ho = 70 kms"^ Mpc"^), Arp 220 belongs 
to a class of galaxies known as Ultraluminous 
IR galaxies (ULIRGs). These types of objects 
are characterized by having an IR luminosity of 



1 



10^2 < L^^^ _ Ljj^ [8 _ 1000 /im] < 10^^ L© . It is 
not entirely clear yet what exactly the source pow- 
ering those high luminosities is. Different sce- 
narios have been considered: The nuclear power 
source may be either a strong nuclear starburst 
( Joseph|[l999[ ), a heavily obscured active galactic 
nucleus (AGN, [Sander s|p!999 ) or a mixture of the 
two. Independent of what the true energy source 
is, a large molecular gas concentration in the cen- 
tral kiloparsecs is needed in order for the energy 



output to be sustained (see e.g. Sanders et al. 



1988). 



Since ULIRGs as a population show a mix- 
ture of starburst and AGN properties, they may 
trace different phases in a ULIRG-to-QSO evolu- 
tionary sequence. Various scenarios linking the 
evolution of ULIRGs to QSOs have been sug- 
gested (e.g. [Sanders et al.|1988[|Bekki et al .|2QQ6D. 
These scenarios allow for a weak, but growing. 



AGN (Taniguchi et al. 1999) at the center of 



Arp 220. Theories competing for an explanation 



of the observed high IR luminosity (Sakamoto 



et al. 2008) and the rather high star formation 
rate of 340 Moyr"^ (|Baan[ |2QQ7[) claim, e.g. a 



pure nuclear starburst contribution, contributions 
from hot cores or a mixture of AGN and star- 
burst. However the detection of Fe Ka emission 
at 6.3 keV at slightly more than 3a significance 
may also indicate the presence of a deeply dust 
embedded AGN. Observations of the [Oiii]5007 
line which is a good AGN tracer but also sensitive 
to star formation, yielded no detection and there- 
fore support the extremely high extinction found 



for Arp 220 (see e.g. [Colina et al.|[2004|) . Opti- 
cal observations reveal that the galaxy is clearly 
crossed by a dust lane spreading from NE to SW 
and show the presence of extended tidal tails ( jArp| 
19661 indicating that the galaxy is undergoing a 
merging event. 

Following it's position in the proposed evolu- 
tionary ULIRG-to-QSO scenarios, Arp 220 can 
provide the link between ULIRGs and elliptical 
galaxies. It also can be considered as a bench- 
mark for understanding the complex properties 
of ULIRGs and the objects representing transi- 
tion stages to QSOs in detail. Arp 220 may also 
turn out to be useful to explain the rotation-like 
structures detected in high-redshift galaxies (see 
Forster Schreiber et al.[|2006|). Studying Arp 220 



resolution configurations may allow us to draw 
conclusions for radio detected galaxies at high 
redshifts which we cannot study in this great de- 
tail, yet. 

High resolution radio and near-infrared imag- 
ing revealed a complex nuclear structure (see e.g., 
Downes fc Eckart[|2007[ ): The core of Arp 220 ap- 



pears to have (at least) two different nuclei (or 
nuclear regions, e.g. |Norris[[T988l ). They are iden- 
tified as Western and Eastern nuclei, separated by 
about 0.4 kpc. CO(l-O) observations revealed an 
underlying rotating kpc - sized molecular gas disk 
( jScoville et al.||1997 ). 

The brightest source in the complex structure 
of Arp 220 in the near-IR ( jScoville et aL][1998D is 
the western nucleus, Arp 2 20- West. The high res- 



olution CO (2-1) map (Scoville et al. 1998) and 



the 1.3 mm continuum maps (Downes & Eckart 
2007[ ) show peaks at the same location. The 



2.2 /im/1.1 /im color map of Scoville et al. (1998) 



however reveals that the peak of the dust distribu- 
tion is not associated with this nucleus ( [Downes fc 



Eckart |[2QQ7| ). [Downes fc Ecka rt' also modeled the 



CO emission by a compact, hotter (170 K) dust 
core surrounded by a cooler (50 K) molecular gas 
ring or disk. They inferred a black body lumi- 
nosity of the compact dust source (~ 10^^ L©), 
which implies that the energetic process powering 
Arp 220- West cannot be of starbursting nature. 

The second brightest source in the complex 
structure of Arp 220 in the near- IR (Scovillel 



et al. 1998) and CO emission maps is the east- 
ern nucleus, Arp 220-East. This nucleus is di- 
vided into two components: North-East (CO-NE) 
and South-East (CO-S E pownes fc Ecka"rtl|2007D . 
Sakamoto et al.[ ( [2009[ ) found a deep absorption in 



the CO-NE/SE region comparable to the absorp- 
tion in the CO emission of Arp 220- West ([Downes 
fc Eckart] |2007[ ). Furthermore, the continuum 



peaks at the position between CO-NE and SE at 
the position of the CO-East peak in the velocity 
dispersion distribution ( Downes fc Eckart] 200 7[ ), 



and other ULIRGs in the local universe in high 



thus implying the region inbetween CO-NE/SE is 
associated with a deeply embedded and/or highly 
dust obscured nucleus of Arp 220-East. 

The velocity fields obtained from the observa- 
tional data of the nuclear region show rotation- 
like patterns at different scales (up to 1 kpc). 
By contrast, observations of the outer regions (at 
scales of ~5 kpc) at optical wavelengths present 
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velocity patterns which are disordered and there- 



fore apparently dominated by the merger (Colina 



et al. 



2004). Influences in the inner regions are 



shown on the 2 kpc scale where the outflows com- 
ing from the dust-enshrouded nucleus affect the 
ionized gas (Arribas et al. 2001). Velocity field 



information determined from observations of the 
CO (1—0) emission line shows a rotational motion 
that corresponds to the overall rotation of the un- 
derlying molecular disk ( Scoville et al.|1997 ). The 
velocity field for CO (2—1) however, looks differ- 
ent: Since the CO (2—1) line transition probes 
higher gas densities, one can look deeper into the 
gaseous environment with this higher molecular 
transition. Therefore, a more complex rotational 
pattern emerges: we find indications for rotation 
in CO- West, CO-NE and SE, with the center lo- 
cated in-between the two intensity peaks (see e.g., 
Genzel et aLl|2QQT| [Mundeh et al.||2QQT| |RBviIqs1 



et al.||2003| |Rodriguez-Rico et al.||2005D . Tak- 
ing the findings for higher resolution CO (2—1) 
data into account, one could also understand the 
larger scale CO (1—0) velocity field as a combi- 
nation of the two rotation patterns of the two 
respective nuclei, but smoothed to a lower res- 
olution. Several models to explain the complex 
dynamics and kinematics in the molecular gas of 
Arp 220 have been presented in the literature e.g.. 



the model of one underlying molecular disk ( Scov- 



ille et al. 1997) with additional gaseous disks in 



the eastern and western nuclei ( [Sakamoto et al. 



1999). Eckart & Downes (2001) tried to explain 



the flux concentrations in Arp 2 20- West and East 
by crowding of inclined circular orbits of a warped 
disk whereas Downes fc Eckart] ( [2007 ) modeled 
their high-resolution CO (2—1) observations for 
Arp 220- West for the presence of a black hole. 
They found indications that the western nucleus, 
in fact, contains a black hole with an accretion 
disk. 

In Sect. 2 we describe the observations and data 
reduction. Sect. 3 reports on the results of the data 
analysis in Arp 220-East. The results of the anal- 
ysis of the large scale structure follow in Sects. 4. 
Sect. 5 concentrates on the dynamical modeling of 
the merger in Arp 220. 

Unless otherwise stated, r^A = 0.7, r^M=0.3, 
and Ho = 70 kms~^ Mpc~^ are assumed through- 
out the paper. 



-0.5 
Offset in RA (arcsec) 



Fig. 1. — East-west (EW) CO (2-1) position- velocity cut 
through Arp 220 from CO-NE to CO-SW, with the 1.3 mm 
continuum subtracted, considering a position between CO- 
SE and CO-NE as the rotation center and a radius of 0.3^^ 



2. Observations and data reduction 

In this paper we work with HST NIR and PdB^] 
CO 1-mm and 3-mm observational data. The line 
and continuum observations of the CO 1 mm and 
3 mm emission lines on different spatial scales 
and resolutions were carried out with the IRAM 
Plateau de Bure interferometer (PdBI) in different 
configurations during the years 1994, 1996, 1997 
and 2006, respectively. The data were reduced 
and analyzed with the help of IRAM's GILDA^ 
software. Depending on the observational setups 
and the data reduction, the synthesized beam sizes 
range from 0.3 0^^ x 0.30^^ to 4.99^^ x 3.50 ^^ The 
CO (1-0) high (Dow nes fc Solomon|1998| and low 
resolution data were taken in winter 1996 with 
a bandwidth of 500 MHz and a channel separa- 
tion of 2.5 MHz, resulting in a velocity range of 
^1300 kms~^ and a velocity channel width of 
7 kms~^. Using natural (high resolution) and uni- 
form weighting (low resolution) in the mapping 
procedure gave beam widths of 1.58^' x 1.11'^ (high 
resolution) and 4.99^' x 3.5^' (low resolution). The 
HST/WFPC2 F814W archive image was reduced 
'on the fly' using the best reference files at the 
moment of retrieval. In this article this image is 
shown only as a reference. Details on the NIC- 
MOS data reduction can be found in I Scoville et al.l 



(1998). 



■"^ Based on observations carried out with the IRAM 
Plateau de Bure Interferometer. IRAM is supported by 
INSU/CNRS (France), MPG (Germany) and ION (Spain). 

2http://www.iram.fr/IRAMFR/GILDAS 
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co-West, 



CO-NE 




Fig. 2. — CO (2—1) emission distribution, velocity field and velocity dispersion distribution of Arp 220. CO (2—1) distribution 
of the central ?>" ^ with a beam of 0.30^^ (left), velocity field (middle) and velocity dispersion distribution (right), showing two 
prominent peaks at the positions of CO-West and CO-East and a possibly fainter third one towards CO-SW. 



3. The kinematics of Arp 220-East 

The eastern nucleus is veiled at optical wave- 
lengths. In the mm wavelengths it can be resolved 



into a NE and a SE part (Downes & Eckart 



Former publications (Downes & Solomon 



Sakamoto et al. 1999), before Downes & Eckart 



2007). 



1998 



(2007), did not resolve the eastern nucleus in two 
components in CO emission. The velocity gradi- 
ent which they used to perform their modeling was 
therefore different from the one shown in Fig.[l]- 
in this paper for the first time the CO intensity 
maps are actually confirmed by the p-v diagram 
(Fig.[l]). In the near-IR, CO (1-0) and CO (2-1) 
emission maps the eastern nucleus of Arp 220 is 
the second brightest source after Arp 220- West. 
Higher resolution CO observations enabled [Scov- 



ille et al.| ( |1997| ) to perform a more detailed com- 
parison between NIR and molecular gas emission: 
The maximum of the near-IR color map (see Fig. 2 
of Scoville et al.|p!998 ) is located at a position be- 
tween CO-SE and CO-NE and hence indicates 
that this position is the peak of the dust distri- 
bution. Using H — K and a standard extinction 
law, the Ay for the SE and NE components in 



the NIR are 24 and 18 mag, respectively ( [Rieke fc 
Lebofsky|[T985| . 



The CO (2-1) spectra for Arp 220-East (Fig.[3| 
show two different velocity components for CO- 
SE and CO-NE. The dominating CO line in 
CO-SE has a zero-velocity of vo~-150 kms~^, 
whereas the CO emission in CO-NE is centered 
at vo'^250 kms~^. The area between CO-SE 
and CO-NE marks a transition region in terms 



of CO line components - both lines are present 
but at lower intensities. The p-v diagram (Fig.[T]) 
for Arp 220-East shows a variation from negative 
('U~-200 kms~^) to positive velocities (v^lOO- 
150 kms~^) considering CO-SE as the center of 
rotation and a radius of about 0.3'^ The CO veloc- 
ity field also shows rotational patterns for the CO- 
SE/CO-NE region and for CO-SW. The velocity 
gradients across the regions are Ai'~200 kms~^ 
(CO-West), Av-400 kms-^ (CO-NE & SE), and 
Av-200 kms-i (CO-SW). The velocity disper- 
sion distribution of Arp 220 (Fig. [2] right) shows 
two very distinct peaks, one at the position of 
CO-West and the other centered in the region 
between CO-SE and CO-NE. The widths of the 
emission lines rise steeply from ~ 200-250 kms~^ 
up to ~400 kms~^ in CO-East and even up to 
^500 kms~^ in CO-West. Taking the contribu- 
tion of the underlying disk (visible in CO (1—0) 
and (2—1) as shown by e.g. Scoville et al.|p!997| 



Downes & Solomon 1998) into account results in 



velocity dispersions well in agreement with values 



from e.g. Scoville et al. (1997). 



With the information provided by the CO spec- 
tra (Fig.[3| the emission line width peak being lo- 
cated in-between CO-SE and CO-NE {Yig^ight) 
and the structure of the velocity field (Fig^ mid- 
dle) it is possible to assert the presence of rotating 
material around the CO-SE/CO-NE region. This 
seems to point out the CO-SE/NE region as the 
choice to be identified with a deeply embedded 
and/or highly dust obscured nucleus of Arp 220- 
East. 

In addition to CO-NE and CO-SE we find that 
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-200 200 



Fig. 3. — CO (2-1) spectra across Arp 220-East. The po- 
sition offsets in RA and DEC are marked for each spectrum. 
The spectra are offset by 0.1^^ from each other. North is 
up and east to the left. The (0,0) position is centered on 
the eastern velocity dispersion peak (Fig.|2|. The veloci- 
ties range from -400 to +400 kms"-*^, the intensity range is 
-0.035 < F < 0.024 Jy/beam. Note the presence of two ve- 
locity components moving from CO-SE (vq ~-150 kms"-*^, 
dashed-dotted line) to CO-NE (vq ~ 250 kms-\ dashed 
line). The bars indicate the widths of the two lines. 



a possible third component, CO-SW (south-west, 
discovered by [Downes fc Eckart 2007), is clearly 
identified in the near-IR fScoville et al.||1998 ) and 
in CO (2—1) (line width enhancement ^ 50 kms~^ 
compared to the surroundings of CO-SW, Fig. [2] 
right) ^ although it is not detected in the 1.3 mm 
continuum (see Fig. 2 of Downes fc Eckart||2007[ ) . 

The extinction values derived for CO-SW from 
the NIR are Ay = 17 mag. The structure of the 
velocity dispersion map (Fig.[2| suggests an addi- 
tional peak towards CO-SW. This enhancement 
could be produced by a mass concentration, or 
alternatively by turbulent motions due to the 
merger process. The velocity field of CO-SW 
shows a rotation pattern centered in this region, 
which is also refiected by the local gradient in 
the p-v diagram (Figs.[l] [2|. Although the true 
nature of CO-SW is still unclear, it could be an 
additional very faint nucleus of a possible mi- 
nor merger component (mergers that involve a 
gas-rich disk galaxy and a bound companion or 
satellite that usually has < 10% of the mass of the 
gas-rich galaxy), a huge star forming region or a 



remnant of the merging process from one of the 
colliding galaxies. 

In the simulations of the merger in Arp 220 (see 
Sect.[5|, the CO-SE and NE peaks are handled as 
the core of merger component II, the western nu- 
cleus is the core of merger component I. 

4. The large scale picture 

The structure of the central region of Arp 220 
has been studied extensively over the years. These 



studies (see e.g., Scoville et al. 


1997 Sakamoto 


et al.||1999t Downes fc Eckart 


2007) have shown 



that the CO and continuum brightness imply a 
total gas mass {H2-\-H) that approaches 10^^ Mq 
and a mean nuclear H2 mass density around sev- 
eral 1000 Mq pc~^. The total mass of the eastern 
and western nucleus is still under discussion. It 
depends critically on whether the dynamical mass 
is totally dominated by the molecular gas or rather 
the stellar and possibly (at least in the western nu- 
cleus) black hole mass contribution. The best one 
can say right now is that the masses of the two 
nuclei (excepting the larger scale molecular disk 
they are located in) may be similar to within a 
factor of two (or more). Little effort, however, has 
been dedicated to study the more external regions 
although there are indications for star formation 
activity in these regions, too. 

The elliptical-like envelope of the underlying 
molecular disk (see e.g. the CO (1—0) velocity 
field map in [Downes fc Solomon 1998) may be 



consistent with a late merger phase within the 
elliptical-through-merger scenario. The f act that 



near-IR light profiles follow an r^/^-law ( Wright 
et al.||1990 ) corroborates this statement. This law 



found by de Vaucouleurs fl953") describes the char- 



acteristic radial dependence of the surface bright- 
ness distribution of elliptical galaxies. 

In this paper we study the molecular gas out- 
side the central 3'^ which has not been subject 
to a detailed analysis in the previous literature. 
In CO (1—0) we find indications for emission ex- 
tending about 10^' to the south as well as ex- 
tended emission to the north and west of the CO-E 
and CO-W components (Fig.|4|. Our recent more 
careful re-analysis of the data used in [Downes fc| 
Solomon (1998) makes us confident, that the de- 



tection is real. A further comparison with other 
low-resolution CO observations yielded the same 
result. 



5 




The molecular mass for the extended emission 
towards the south of Arp 220 is determined from: 



15 34 58.0 57:5 



Fig. 4. — Integrated intensity map of the CO (1-0) 
emission in Arp 220 integrated over a velocity range of 
730 kms"-*^ with a beam size of 1.57^^ x l.ll^^and with con- 
tours starting at the la level with increasing intensity in 
steps of 2a. 




Fig. 5. — CO (1—0) emission distribution smoothed con- 
volved with a 7.8^^ beam, centered on the position of 
Arp 220. Positive and negative contour levels (indicated 
in red and blue) at steps of 1.25cr, excluding 0. The grey- 
scale indicates the intensity. The bar represents a size of 
roughly 10 kpc. 



M ={N 



H2 
X 



)-A-2-ii-mn Mq 



(1) 



M is the total mass of molecular hydrogen of 
the extended emission region, {N{H2)) is the hy- 
drogen column density determined from the low- 
resolution CO (1—0) data, A is the area projected 
onto the sky in the distance Di, in cm~^, mn is 
the mass of one hydrogen atom (1. 67x10"^^ kg) 
and /i is the average molecular weight per hydro- 
gen atom (/i = 1.36). Applying this formula to the 
CO data for this region results in a total molecu- 
lar gas mass (within the la level contour line) of 
-1 X 10^ Mq. 

In Fig. [5] we show the line intensity map con- 
volved with a 7.8^' circular beam. The southern 
component appears as a natural extension of the 
central source. A few shallow negative regions sur- 
round the source plus southern extension symmet- 
rically at a radius of about 20^' followed by a pos- 
itive noise belt at a radius of 30" . 

The main source plus its southern extension lies 
at the center of those rings. The symmetry of the 
first negative ring supports the significance of the 
extended positive source structure in the center. 
At that resolution, at which the source size of the 
southern extension is well met, the positive south- 
ern extension reaches 5<j. 

5. Dynamical modeling of the merger 

Arp 220 is the archetypical infrared lumi- 
nous galaxy merger. It shows tidal tails and 
a disturbed central region hosting two promi- 
nent nuclear sources. Motivated by the detec- 
tion of molecular gas in the more extended host- 
galaxy regions of Arp 220 and the two nuclei, we 
present a merger model based on the collision- 
less, non-selfgravitating multiparticle code Iden- 
tikit ( Barnes fc Hibbard||2009 ) . Before describing 
the details of the modeling we justify why a simu- 
lation with a collisionless multi- particle code can 
be used despite the apparently large mass of dense 
molecular material at the merger nucleus. 

5.1. Motivation 

While the nuclear region of Arp 220 contains 
significant amounts of dense molecular material, it 
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encompasses only a few per mill at most of the en- 
tire volume over which the merger takes place and 
has been modeled by us. For the dominant por- 
tion of the host we assume as a working hypothe- 
sis that across the Arp 220 merger the molecular 
gas predominantly occurs in clumps with sizes of 
about 1 pc (sizes of 10^^^ pc are quoted by e.g. 



Parmentier et al.|2Qll||Fontani et al.|2QQ2|[HQEer 



et al.||2QQQ||Vallee ||2QQQD T 

With the help of the mass of the gas clumps 
and the total gas mass of the extended region to 
the south of Arp 220 the particle collision rate can 
be estimated. Using the critical density commonly 
found in gas clumps (ncnt = 10^ ...10^ cm~^), the 
clump mass (Mdump) can be derived by 



Mel 



ump 



(2) 



V is the volume of one clump (~ 1 pc^), mn 
(1.67 X 10~^^ kg) is the mass of one H atom and 
^ciump is given in units of solar masses Mq. Since 
we calculated the gas mass of the southern re- 
gion of extended CO emission (Mgas ~ 10^ M©), 
the clump density can now be determined, via 



Mclump / NdV = M, 



gas 



(3) 



where Mdump and Mgas are known, N is the clump 
density and dV is the volume of the southern re- 
gion of extended CO emission. With the ve- 
locity dispersion a ('^10 kms~^) and the cross- 
section s (^ 1 pc^) known, finally the collision rate 
i^coii can be derived by: 



^coll = N G S. 



(4) 



From this value the time tcoii between two colli- 
sions is then determined by the reciprocal of the 
collision rate. This value gives an upper limit for 
tcoii of particles. Values for the time scale range 
between 3.6 x 10^ and 3.6 x 10^^ yrs, dependent 
on the critical density of the gas clumps used in 
the calculations. These time scales are compara- 
ble to the dynamic time scale of the merger itself. 
The other values for the time between collisions 
are even larger. Even if agglomeration between 
the gas clumps should occur, the energy contribu- 
tion from the merger process, via tidal forces, to 
the velocity dispersion of the gas is large enough 
that the clumps in the outer regions do not col- 
lide any more. Therefore the gas clumps can be 



treated, to first order, as stars. Both, stars and 
gas clumps can then be handled, to first order, as 
collisionless particles. This fact is implemented in 
the Identikit code, so that it is justifiable to use 
this modeling tool for the purposes of the simu- 
lation presented. For that particular reason the 
Identikit model can be used without restriction. 

Also the contribution of the atomic gas compo- 
nent can be neglected: Since Arp 220 is already 
in the final stage of merger the interacting galaxy 
component particles already collided, in the cen- 
ter of Arp 220, at least 1-2 times with each other. 
Therefore the HI is either stripped from the region 
studied with the merger simulations (Fig. 3 of Hi- 



bbard et al. 2000) or is already transformed into 



molecular gas (H2). 



5.2. The Identikit model 



The Identikit model ( [Barnes fc Hibbard|[2QQ9l ) 
is a tool to simulate mergers of galactic disks with 
the help of test particles. Test particles (represent- 
ing stars and clumps in this simulations) are used 
due to their ability to reproduce features such as 
bridges, tails and shells, in particular, very well. 
The simulations described by the Identikit model 
are based on collisionless N-body simulations (for 
the validity of the assumption of collisionless parti- 



cles see Sect. 5.4). The simulations presented here 



are based on a model where a bulge, a disk, and a 
halo component are included. 

The initial orientation of the two galaxy disks 
related to the orbit is described by the inclination 
i and the peri-centric argument uj of each galaxy 
disk. Parameters defining the orbit are the eccen- 
tricity of the orbit e, the peri-centric separation p 
and the mass ratio /i of the two disks. In the Iden- 
tikit model the eccentricity is fixed at a value of 
one. The evolution time t describes the age of the 
merger in internal units of— 2<t<8. ^ = repre- 
sents the point in the merger evolution at which, in 
the configuration of the merging/interacting sys- 
tem, the idealized Keplerian orbit reaches the peri- 
center. 
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a) c) 




Fig. 6. — Identikit results for the model with the best fit to the available data of Arp 220. a) to d) show the model projected 
onto the x-y (a), x-v (b), v-y (c) and x-z (d) planes. In e) the magnified x-y plane is shown. The grey bar represents a size of 
roughly 10 kpc and the black bar represents a velocity scale of ~ 500 kms"-*^. 



5.3. Simulation results 



Optical HST archive composite imag€ ^ and 



mm- CO (obtained by Eckart fc Downes| |2QQ1 
Scoville et al.|p!997 ) data are used to fit the Iden 
tikit model simulations. The HST images were 
used to identify the 'overall shape' (x-y plane) and 
the CO (2—1) data were used to determine a mea- 
sure on the position- velocity diagrams (x-v plane). 

In order to obtain a grid of possible models and 
to obtain the 'best model fit' from that grid, many 
different simulations were run systematically. To 
get a first estimate on the model parameters for 
Arp 220, the model of the Antennae Galaxy was 
used as a starting point, especially for the in- 
clinations i and pericentric arguments uj of the 
galaxy disks (Antennae: ii=i2= 60°, ui=U2=- 
30°,| Toomre fc Toomre||1972[ [Barnes||1988| ). The 



Antennae Galaxy has a similar morphology, com- 
pared to Arp 220, and it also shows tidal tails from 
both merging galaxies, but the viewing angle is dif- 
ferent. The eccentricity was fixed internally from 
the start at a value of e = 1. In the Identikit sim- 
ulation tool, only two values for the mass ratio /i 
were available to choose from, 1:1 and 1:2 with 



^from http://imgsrc.huhhlesite.org/hu/dh/images/hs- 
2008-16-aq-fulLjpg.jpg (Credit: NASA, ESA, the 
Huhhle Heritage (STScI/AURA)-ESA/Huhhle col- 
lahoration, and A. Evans (University of Virginia, 
Charlottesville/ N RAO /Stony Brook University)) 



the eastern nucleus being the mass component 1. 
Actually, the simulations were performed for both 
mass ratios. At the start of the simulations all 
possible values for the pericentric separation p, in 
the combination with the merger time scale t, were 
run. After the first few crudely determined simu- 
lations, the parameters, e, /i, p, ii, c^i, ^2, ^2 and 

were better constrained to a smaller parame- 
ter space, the viewing angles were varied. First 
only one of the angles was changed, the other two 
remained fixed, later all were changed in combi- 
nation. The variation of the viewing angles was 
performed systematically in steps of 10°. 

Changing the parameters and the combinations 
thereof resulted in a set of parameters describing 
the model that fits the observational data best 
(Fig.|6]). This best fitting model was obtained by 
a visual inspection of the model fits and the ob- 
servational data and by applying a set of different 
criteria (see also Fig.[9|. This best fitting model 
has a point in the time evolution that results in a 
good match between the position and velocity of 
the 2 nuclei (Figs.[6j [t] and |8| and represents es- 
sential features of the more extended distribution 
of the molecular gas ( Figs. [t] and |9|. 
The set of Identikit model parameters delivering 
the best solution fitting of the HST image and the 
CO position-velocity diagram is comprised of the 
following parameters: The eccentricity at a value 
of e = 1. For the pericentric separation p a value 
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Fig. 7. — Arp 220 HST image overlay with the best model. Left: HST ACS/WFC F435W and F814W composite image 
of Arp 220 with an increased contrast in order to show the tidal tails of the galaxy more clearly. Right: The same, but 
superimposed with the x-y plane of the best model fit and also with an increased contrast (between model and background 
image). The white bar represents a size of roughly 10 kpc. 



of 0.125 delivered the best results. A mass ratio 
of the two galaxy disks of 1:2, an unequal merger 
(see Figs. [g] and [t]) resulted in the best model fit. 
The lower mass galaxy (see Fig. [g] in magenta) was 
initiated with an inclination i of 60° and a pericen- 
tric argument uj of —15°, the higher mass galaxy 
was described initially by an inclination i of —60° 
and a pericentric argument cu of 40° . The viewing 
angles of Ox = 90°, Oy = 180° and 6z = 40° gave 
a good match of the simulation data to the obser- 
vations. The evolution time scale for the merger 
resulted in a value of t = 3 (in internal units), la 
errors of the Identikit modeling parameters are of 
the order of 8.3° for the viewing angles 9x, Oy 
and The errors for the internal time scale, 

the inclinations and pericentric arguments of the 
disks are At^O.l, Ai^ 10° and Auj --10°. A com- 
parison of the obtained model fit of Arp 220 to 



the study of Gonzalez- Garcia & Balcells (2005), 



of merging galaxy systems, results in a transla- 
tion of the evolution time t to a merger age of 
roughly ~ 6 x 10^ yrs. This value is in very good 
agreement with the age of the merging system de- 



termined by Mundell et al. (2001). They give a 
value of 7 X 10^ yrs. 

The comparison of the HST image and the best 
model fit in Fig.[7| shows that the model very well 



represents the properties of the Arp 220 in the 
HST observations. The sharp edge towards the 
north-east is well represented by the model. The 
tidal tail of the less massive galaxy (in magenta) is 
in agreement with the extended emission towards 
the north-west. The more southern tidal tail is, at 
the basis towards the central region of Arp 220, in 
good agreement with the HST image. Further out 
this extended emission does not coincide with the 
model any more. As shown in Fig.jSj the overall 
shape of the model matches the CO p-v diagram. 
Note that a difference is that the CO p-v diagram 
was taken in the central region of Arp 220. Hence 
it does not cover the same velocity range as the 
simulations do and does not take the regions fur- 
thest outwards into account (see Fig.[8|. The ex- 
tensions to the NE and the SW side of the zero- 
offset axis of the p-v diagram describe the CO p-v 
diagram very well. The difference in the veloc- 
ity range originates from the test particles further 
outward of the central region not covered by the 
CO data. Also the extended emission found in the 
low resolution CO observations (Fig.|4| coincides 
very well with the simulations: At the location of 
the extended emission ~ 10'^ to the south of the 
central part of Arp 220 a high test particle den- 
sity is found in the Identikit simulations (Fig.[6|. 
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Fig. 8. — Comparison of the posit ion- velocity diagram obtained from the Identikit simulations (left) with the p-v diagram 
from CO (2—1) observations (right). The quality test parameters n (in green), e (in dark blue) and 6 (in red) are indicated in 
both plots. Note that the CO p-v diagram does not cover the same velocity range as the simulations do, since it was obtained 
for the center of Arp 220 and does not take the regions furthest outwards into account. 



It is also worth noting that in Fig. [S] the density of 
the simulation dots on the left does neither reflect 
the molecular gas mass solely nor does it reflect 
the emissivity of the molecular gas. Also the right 
hand image is spatially convolved with the beam, 
whereas the left hand site is not convolved to allow 
for a clearer view onto the distribution of simula- 
tion dots. 

The quality of the model fit was determined 
by the following criteria: Size and compactness of 
the merger galaxy, and by the location and posi- 
tion angle of the tidal tails. The position angles, 
a and /3, of the tidal tails to the north-west (rep- 
resented by the yellow line in Fig.[9| and closer to 
the central region (represented by the green line 
in the same Figure) of Arp 220 have values of 
a - (124.6 ± 9.6)° and /3 - (6.9 ± 7.6)°. The sharp 
edge of the central region of Arp 220 is located at 
a position angle 7 of ~ (125.2 ± 4.8)°. The surface 
filling factor, a representation of the size and com- 
pactness of the galactic central bulge, was found 
to be F-(79±6)%. This filling factor F rep- 
resents the area within the box, shown in Fig.[9j 
occupied by the central region of Arp 220. F is 
determined by fraction of the area of the box and 



the box area, with the triangles subtracted where 
the particle density is reduced to about zero: 



F = 



box area — triangles 
box area 



% 



(5) 



All errors stated for the results are errors at the 
la level. By direct comparison, this la error for 
each criteria parameter was determined from the 
HST image. For the error values resulting thereof 
a (Eq. 3.1) of 1.3 was derived by 



A^obs 

E 

i=l 



(/o 



bs,i 



/model, i) 



(6) 



as a measure for the deviation of the best model 
fit from the observed data. /obs,i is the respective 
criteria parameter measured in the observational 
data, /model, i IS the parameter determined from 
the best model fit and a^ is the observational er- 
ror. Assuming Ax^ ~ 1 for the variation of each 
Identikit parameter separately (one parameter was 
variable, the others were kept fixed), the internal 
errors for the Identikit model parameters were de- 
rived (see before). 10 model parameters in the 
Identikit tool are faced by 10 criteria parameters 
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Fig. 9. — The model evaluation criteria are shown on top of the HST image (left) and the best model fit (right) for a 
comparison of the observed data to the model fit and an assessment of the quality of this model, a, (3 and 7 represent the 
position angles of the two tidal tails and the sharp edge of Arp 220 to the east, respectively. F is the surface filling factor, 
which is a measure for the size and compactness of the galaxy. The bar represents a size of roughly 10 kpc. 



obtained from the observations: location (six pa- 
rameters) and position angle (three parameters) of 
the tidal tails and the sharp edge of Arp 220, and 
the surface filling factor (one parameter). With 
the help of these criteria the quality of the fitted 
models was assessed. The comparison of the best 
model to the HST image in Fig. [9] shows how well 
it fits the observational data. 

We furthermore check the quality of the iden- 
tikit model fit versus the CO p-v diagram. We 
therefore use three parameters to compare the CO 
p-v diagram to the p-v diagram obtained from the 
identikit simulations: n (in green) represents the 
distance between the two nuclei, e (dark blue) is 
a measure for the distance between the NE and 
SW extensions, relative to the zero-offset axis and 
5 (red) describes the angle between the zero-offset 
axis and the axis on which the two nuclei are situ- 
ated. All three parameters are in good agreement 
within the error bars, defined by the error of the 
nuclear position in the model and the scatter of 
the off-axis points in the two extensions with re- 
spect to the zero-offset axis. 

5.4. Discussion 

Since it is known from observations at mm 
wavelengths that the mass ratio of the two nu- 
clei Arp 220- West and Arp 220-East is about 1:1 
(see e.g., Scoville et al.|1997 Sakamoto et al.|1999 



Downes fc Eckart|2007D , it is curious that the sim- 
ulation results presented here give a value for the 
mass ratio fi of the parent galaxies of 1:2. One pos- 
sible explanation is that the model, is not unique, 
i.e. the set of parameters fitting the HST and 
CO data. Nonetheless, the best model obtained 
here produces a 'fair' fit to the HST image and 
the CO data. Another point that speaks in the fa- 
vor, of the model presented here, is the determined 
merger age that is in very good agreement with the 



literature value from Mundell et al. (2001). 



However, the situation may be more compli- 
cated since (and this can be regarded as a short 
coming of our fitting procedure) the data that were 
used to fit the model to, were taken at two differ- 
ent wavelengths and are therefore basically tracing 
matter under different physical conditions. Only 
in the central arcseconds the gas has an impact on 
the overall gravitational potential, comparable to 
the stars impact on the potential. Further out the 
stellar potential clearly is dominating. Therefore 
the gas follows the potential of the stars and con- 
sequently the stars themselves. Since it was easier 
to guide the fits from the simulations by the tidal 
tails, seen in the HST image, more weight is put 
on the outer regions of the galaxy. Hence it is 
justifiable to use the HST data as well as the CO 
observations to fit the model results to. Outside 
the central region of Arp 220 gas is assumed to be 
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counterexample 1 



counterexample 2 



observational data 





Fig. 10. — Counterexamples to the best model fit. Identikit simulation with the same parameter set as in the best model fit, 
but with a mass ratio fi of 1:1 (left) and a viewing angle Bx of 0° (middle). The x-y plane is shown at the top (a) and c)), the 
positi on-velocity diagram (x-v ) at the bottom (b)and d)). The HST composite image and the CO position- velocity diagram 
(from [Eckart Sz Downes|2001| used to fit the Identikit simulations are shown in e) and f) for comparison. Note that the HST 
image is rotated by 180° to match the Identikit internal axis orientation. The grey bar represents a size of roughly 10 kpc and 
the black bar represents a velocity scale of ~ 500 kms~^. 



found in clumps with sizes of about 1 pc. There- 
fore the gas clumps can be treated, to first order, 
as stars. Both, stars and gas clumps can then be 
handled, to first order, as collisionless particles. 

To show that the model based on the dynamical 
multi-particle simulation presented here is indeed 
one possible qualitative description of the system, 
two possible counterexamples are shown for com- 
parison in Fig.[lOj Since a mass ratio of 1:1 be- 
tween the nuclei is commonly assumed, one coun- 
terexample representing this mass ratio is shown 



in Fig. 10 (left). It was modeled with the same set 
of parameters as for the best model fit, but with 
the different mass ratio ja. In this figure it is very 
clear that the model does not fit either the HST 
image nor the CO p-v diagram (Fig.[lO|. The same 
applies to the second counterexample (Fig.jlOj in 
the middle). Here the same set of parameters as 
for the best model fit was used again, only this 
time with a viewing angle Ox of 0°, which cor- 
responds to an angular deviation from the best 
fitting model of lOa. The resulting model very 



clearly does not fit the underlying observational 
data. 

6. Conclusions 

In this paper on Arp 220 the results on the 
study of the interferometric CO data and the Iden- 
tikit simulations of the merger are discussed. The 
eastern nucleus of Arp 220 is studied and the ob- 
ject as a whole is studied for more extended struc- 
ture in CO gas emission. The analysis of the CO 
lower resolution data shows that the velocity dis- 
persion component centered between CO-SE and 
CO-NE is a very promising candidate to be the 
'real', but deeply embedded/highly dust obscured, 
nucleus of Arp 220-East. Indications for emission 
~ 10^' towards the south, as well as to the north 
and to the west of the two nuclei were found in 
the low resolution CO (1—0) maps. Furthermore 
simulations of the merger in Arp 220 were per- 
formed with the Identikit modeling tool. The 
model parameters describing the galaxy merger 
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best give a mass ratio of 1:2 and result in a merger 
of ^6x10^ yrs, which is in good agreement with 
values from the literature, such as Mundell et al.l 



(2001). Apparently Arp 220 can be reckoned a 



very interesting case of merger galaxy, where the 
object can be placed between a minor merger 
(from the CO observations) in one extreme and 
an equal (1:1) merger (from the simulations, CO 
data and published models in the literature) in the 
other. 
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